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I. INTRODUCTION 

Earlier this decade a "low-density theorem" (LDT) for the linear term in the baryon den- 
sity dependence of the quark condensate was formulated |l|-^ . This first term in the density 
power series expansion of the quark condensate is proportional to the nucleon/constituent 
quark S term, which leads one to believe that, perhaps, it is only the chiral symmetry 
breaking of the (strong) hadronic interactions by the current quark masses that controls the 
baryon density dependence of the quark condensate and thus makes it essentially unique. 
Soon thereafter it became clear, however, that the extrapolation of this linear formula from 
the low-density region upwards, to the neutron star densities, would be highly contentious. 

In this paper: (1) We present an extension of the Drukarev-Levin low-density theorem 
to terms of arbitrarily high order in the McLaurin expansion of the fermion condensate's 
dependence on the baryon density p. The n-th McLaurin coefficient is just the connected n- 
nucleon S term (elastic matrix element), which is a function of the explicit chiral symmetry 
breaking (xSB) terms in the Hamiltonian of the theory. (2) We show how the "theorem" 
works in an explicit example: We calculate the 0{p^) coefficient in the lowest perturbative 
approximation to the linear sigma model and then iterate these diagrams to infinite order. 
The result can be resummed and we obtain a closed-form solution - we use this result to 
establish the range of validity of such a diagrammatic calculational scheme. 

II. BARYON DENSITY DEPENDENCE OF THE S TERM 

A. Proposition 

We propose an extension of the low-density theorem of Drukarev and Levin [^ to arbi- 
trary (positive integer) powers of the baryon number/quark density p: 

The n^^ coefficients in the McLaurin expansion of (S)p is the n-fermion/nucleon/quark 
matrix element of the (pion) S double commutator 



-J a=l 

is 

(p |S| p) = 5] ^{jN |S| jiV)eo„nected , (2) 

where the sum over j involves connected matrix elements. 

As preconditions we assume that: (i) the color symmetry of the theory is neither sponta- 
neously nor explicitly broken; (ii) the theory can be described by a Lorentz invariant, local, 
approximately chiral-invariant Lagrangian density with a partially conserved axial Nother 
current; (iii) the axial anomaly does not invalidate the relevant axial Ward identities. 

The proof is based on the current-algebraic relation/identity: 

{T)^ = {a\i:\a) = W{a\[QUQinMm\<^) 

-^ a=l 

= -\fl\iuvJ2{n\k)a\S\7,'^{k)a) , (3) 

between the exact elastic soft-pion vra scattering (S) matrix and the corresponding Heisen- 
berg representation (pion) sigma term matrix element. Here a is the flavour index of the 
external pions, which is averaged over the three varieties of pions, i.e., a = 1,2,3. This 
result can be adapted to the kaon sigma term by specifying a = 4, 5, 6, or 7, or by averaging 
over some subset thereof. Formula (|^) is derived from a chiral Ward identity and the LSZ 
reduction formula for the case of a nucleon (a = A^) on p. 131 - 7 of Ref. 0. [This is 
equivalent to applying Sakurai's "master formula", p. Ill - 2 in Ref. 0], twice, once to the 
initial and once to the final state.] Note that the sigma term is only sensitive to the chiral 
symmetry breaking (xSB) terms in the Hamiltonian density H = H^ + Ti-^sB 

[Qt,[Qt,n]] = [Qt,[Qin^sB]], (4) 

the chiral charges being constants of the motion in the chiral limit 

[Qt,n^]=o. (5) 



Now apply this equation to the nuclear matter state ket |a) = \p) = \nN). We shall 
work in a large but finite box of volume Q in order to avoid momentum-conserving Dirac 
delta functions, which are replaced by Kronecker ones. This also means that we can have 
a finite baryon density p = f| without taking the number of baryons to infinity, n ^ oo, 
which we leave for the very last step. Therefore 

(S), = (p|S|p) 

= lim (nA^|S|nA^) 

= hm \J2{nN\[Ql[QinMm\nN) 

1 ^ 

= — /^ lim limy2{Tr''{k)nN\S\n''{k)nN) , (6) 

a=l 

where we used the fact that the finite density (zero temperature) ground state wave function 
is just the wave function of n fermions/nucleons at rest enclosed in a box of volume Q. 
Formula (^) relates the finite density S term matrix element to the exact (forward) scattering 
amplitude of a soft, i.e., vanishing four-momentum {k = 0) pion from n nucleons at rest. 

Here we have merely used the definition of the S-term to rewrite the object of interest 
in terms of an S-matrix element. To be sure, no new information was gained in this step 
and at no loss of generality. The result is that we may now use the general results of S- 
matrix theory. That, in turn, allows a Feynman- diagrammatic organization of our problem. 
The most important property of the S-matrix in this regard is its decomposability into 
disconnected and connected parts. The latter part falls further into reducible and one-, 
two-, three-particle, etc. irreducible classes. Specifically, in our case of pions and nucleons 
we may separate the exact (Heisenberg representation) S-matrix into the following distinct 
categories: (i) the completely disconnected graph (one pion and n baryons all propagating 
without interaction); (ii) n simply connected vrA^ scattering amplitudes multiplied by {n— 1) 



disconnected baryon lines; (iii) 
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n(n— 1) 



simply connected ttNN scattering amplitudes 



multiphed by {n—2) disconnected baryon line, etc. 0. Each of the disconnected subdiagrams 



has its own momentuiii- conserving delta function, which translates into one volume (fl) 
factor for each disconnected baryon line [] [0. Therefore 
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Now turn this back into a statement about the E terms 

{nN\L\nN) 



{nN\nN) 
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— ) (jiV|S|jW) connected • 



(8) 



We are now ready to take the thermodynamic limit fi ^ oo at (finite) constant density p. 
Since we assumed a large box to begin with, the baryon number n = pQ contained therein 
may be taken to be far larger than any given order in the power expansion, n ^ j, and in 
the thermodynamic limit n -^ oo one may write 

fA 
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n\ 



n-' 



which leads to the announced result 
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(10) 



^This, of course, means that the dimension of the various S term matrix elements varies with the 
number of nucleons (baryons/quarks more generally): the vacuum S term has dimension M^, the 
(single-)nucleon S term has dimension M^, etc., i.e., M^~^" for the n-nucleon S term 



when the nuclear ground state is normahzed (p|p) = 1. It is manifest that this result does 
not hold in finite systems at orders j comparable to the total baryon number n. 

A reminder seems in place that this series need not uniquely determine (S)p for arbitrary 
values of p because its radius of convergence may be small, or even zero. One possible reason 
for such a behaviour of a power series is that it really is a Laurent one, i.e., the function may 
have poles and/or branch cuts in the complex baryon density (p) plane. Explicit calculations 
will show this to be the case. 

B. Consequences 

Note that the "mechanical" sigma operator 

^ ^ a a 

Smech = o 2^ [Qs) [Qs; '^mcch(O)]] = 7imcch(0), (11) 

-J a=l 

equals the "mechanical" xSB Hamiltonian 

'^^mech = ^m^'^ = m^uu + ui'^dd, 

due to the current quark masses, or the bare nucleon mass term 

n^ech = NM^^N = M^pp + Minn. 

Our result Eq. (^) holds either for quarks, or for baryons as the fermions. Since the chiral 
symmetry breaking is defined in terms of quark degrees of freedom, but we do not know 
the solution to the quark dynamics in baryons, we must first establish, or rather postulate 
a relation between the quarks' and the observable hadrons' properties. One assumption 
common in the literature is that Sat = Sag. Assuming, by the same token, effective 
proportionality of the nucleon and quark condensates, our Eq. (|1^) leads to a statement 
about the quark condensate at finite baryon density 



^ It is clear that this amounts to the "impulse approximation", but it was recently shown that 
the two-quark operator corrections must exist M] if chiral symmetry is to be preserved. 



TTOT = P^ = l- jj^, E ^{JN |S| jiV)„.,., . (12) 

(A^A^)o {qq)o [Um^r J^i j' 

The j = 1 term together with the unity on the right-hand side form the Drukarev-Levin 
LDT. The rest of this formula appears to be new 0. 

This form of the LDT has been used to argue about the behavior of the p.s. meson 
masses at non-zero baryon density. We remind the reader that the quark condensate is not 
the only density-dependent variable influencing the p.s. meson mass - one must take into 
account the (model dependent) p.s. meson decay constant's density dependence as well, 
which is not commonly done. Lest this result leads one to believe that the problem of the 
density dependence of p.s. meson masses is solved, we remind the reader that neither the 
radius of convergence of the series is known, nor do we have, as yet, an efficient algorithm 
for calculating the higher-order coefficients in the expansion. 

We next present an example of a perturbative calculation in the (chiral) linear sigma 
model and its iteration to arbitrarily high order. 

III. AN EXAMPLE 

Before we proceed we remind the reader of the assumptions of the LDT: the matrix 
elements must be either the exact Heisenberg representation ones, or approximate matrix 
elements that satisfy the underlying chiral Ward identity (master formula). The former, is 
beyond our powers, the latter are available, mostly in the form of perturbative solutions, 
though solutions to several models have been found that sum certain infinite classes of 
Feynman diagrams (we shall not consider those approximations here). We shall calculate 
the terms of 0{p'^) in the chiral linear sigma model, the 0{p) calculations having been done 
before 0-^. This calculation will be with (free) nucleon degrees of freedoms, new methods 



^ The idea that the elastic vr — nN scattering amphtude is related to the p-dependence of {qq)p is 



implicit in Refs. [§JT§. 



for calculation of sigma terms in interacting nucleon systems having been developed only 
recently [0]. 

A. One- and two-nucleon sigma terms in the linear sigma model 

The one-nucleon connected sigma term can be extracted from the elastic ttN scattering 
amplitude, which in the linear a model is given by the sum of the three diagrams in Fig. 
1: two nucleon-pole diagrams, (b) and (c), and one a-exchange graph (a). These have been 



calculated many times II^] and will not be repeated here. 

The pion-two-nucleon connected matrix element is substantially more complex: there are 
more than 50 connected diagrams in the lowest-order perturbative (Born) approximation. 
All of them are, by definition, reducible diagrams. We separate these graphs into three 
distinct subsets: 1) pion "rescattering" ; 2) pion scattering on a sigma-meson-in-fiight; 3) 
initial-, final-, and intermediate-state interaction diagrams. 

1) Eighteen (2 x 3 x 3) of these graphs form the simplest connected, reducible "pion- 
rescattering" amplitude that is built up from two pion-one- nucleon effective vertices. Fig. 1. 
Their contribution equals 

(.«.^-2(f^)\ (13, 

The currently accepted value of the nucleon sigma term (S)7v lies between 45 and 65 MeV. 

2) Only the sigma terms of scalar states, discrete or continuum, exchanged between the 
two nucleons contribute to the two-nucleon sigma term. [This follows from the fact that the 
sigma term is a Lorentz scalar.] In the first Born approximation to the linear sigma model 
the only such object is the sigma meson, so we end up having to calculate its sigma term 
and its contribution to the the two-nucleon sigma term. Four Feynman diagrams, see Fig. 
2, contribute to the elastic pion-sigma meson scattering amplitude and hence to also to its 
sigma term. The result is 

((T|S|a) = 3m2 . (14) 



Inserting this result into the two-nucleon signia term we find 



2 ' 



2, m 

mi 



(E)^^ = 3<7oM -J ' (15) 



where go = {qaM/ f^,) = 12.6 is the a-nucleon couphng constant. The currently accepted 



value of the a mass lies between 400 MeV and 1200 MeV |T2 . 

3) These are not all of the reducible graphs at this order, however: all of the remaining 
graphs can be described as either initial-, intermediate-, or final-state interactions, see Fig. 
3. They are formed either by having a complete pion-nucleon scattering amplitude attached 
to one external nucleon line, or by having a pion-nucleon vertex attached to two external 
nucleons in the sum of one meson exchanges between the two nucleons (the NN potential). 
This one-meson exchange potential ought to be summed into an infinite ladder by iteration 
of the Bethe-Salpeter equation and thus made to form a proper interacting NN (scattering or 
bound) state. In the end, this summation is equivalent to the introduction of NN correlations 
into the nuclear matter. These diagrams diverge in vacuo, although they are Born diagrams, 
because the intermediate nucleons are on their mass shells, due to the external soft-pion 
condition. 

Points 1) and 2) put together lead to 

(S)2Ar = (S)2Af + {^)2N 

= -2f^)%3f»2!^)\ (16) 

for the total. With the currently accepted value of the a mass, the second term exceeds the 
first one. 

One feature of this result stands out: the appearance of the second power of the single- 
nucleon sigma term coming from the repetition of the one-pion reducible diagrams involving 
the "primitive" (elementary) single-nucleon sigma term contributions. It has been pointed 
out by M. Ericson p[ that this sort of behaviour generalizes to higher orders of perturbation 
theory. Fig. 4, and that the graphs can be resummed into a geometric progression. We shall 
address that question next. 



B. Iteration of primitive contributions in the N-nucleon sigma term 

In the foregone analysis one will have noticed that certain lower-order graphs are repeated 
in the higher ones with simple regularity. The multiplicities of these graphs can be figured out 
and the infinite (sub-)series resummed. We shall show that several classes of such diagrams, 
sum into geometric series. That result is particularly useful because its properties, such as 
the radius of convergence and the analytic structure are well known. Yet, as we shall show, 
uncritical application may sometimes lead to erroneous conclusions. 

Single-nucleon pion rescattering The connected but reducible diagrams. Fig. 4, consist- 
ing of individual vrA^ elastic scattering amplitudes connected by a single pion line carrying 
zero four-momentum can be written as 

{nN |S| niV)rescatt. = n!(S)^ [-mlflf''"'^ , (17) 

where the factor -^^ comes from the propagation of a soft tt, the (S)^?/,^^ comes from 
a single irreducible elastic nN scattering amplitude (vertex) for soft pions, and n\ is the 
number of identical graphs. Hence the infinite series of such connected one-pion-reducible 
graphs is readily summed up as 



{NN), 



{(l(l)(- 



E 



-p{^)n 



n 



(18) 



^^ p{^)n 



if7Tm^y_ 

Since this is a geometric progression, the series has a finite radius of convergence (in p) given 
by 

, I I (/tt^^tt)^ 

P< \pci\ ■ 



{^)n ' 
which equals roughly twice (1.8), or three times the normal nuclear density for (S)jv = 65, 

45 MeV, respectively. 

It seems obvious that the complete series including all diagrams, at best can have the 

convergence properties of the worst-behaved sub-series. In other words, our resummation of 

the one-nucleon pion rescattering diagrams seems to tell us the maximum reliable density 
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calculable with such diagrammatic methods. An important practical consequence of the 
above reasoning and of the empirical value of the kaon-nucleon sigma term (S)xAr = 200 
MeV is that the radius of convergence of the series relevant to the K -condensation is at best 
a fraction of one normal nuclear density. Surprisingly, this conclusion turns out to be only 
qualitatively correct, due to cancellations from other two-nucleon sigma term contributions. 
We shall explicitly show that in the linear sigma model the inclusion of another two-nucleon 
sigma term into the infinite series leads to an increase in the number of poles, change of 
their positions and to an increase in the radius of convergence of the series. 

Two-nucleon pion rescattering The connected two-nucleon "pion-rescattering" graphs 
fall into one of the two categories: (i) even-n, and (ii) odd-n, and will be dealt with accord- 
ingly. 

(i) In the former case, one part of the amplitude is in the form of a product of n/2 
pion-two-nucleon graphs connected by n/2 pion propagators. Since one can select the first 



pair m 
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\^i 



many ways, the second pair in 



^„-2^ 



many ways, etc., one finds that 
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V 2m2/2 , 

The iteration of such two-nucleon pion rescattering diagrams, together with the vacuum 
term, leads once again to a geometric progression: 
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(20) 



2if^m^)\ 

(ii) For n-odd the amplitude is in the form of a product of one vrA^, [or irmN graph, with 
m-odd] and {n — l)/2, [or {n — m)/2] two-nucleon irreducible graphs connected by (n+ l)/2 



11 



pion propagators. Since one can select the single vrA^ graph to be any of the n nucleon lines, 

( n-l \ 



and the first pair in 
one finds that 

(nA^|S|nA^)| 



' n-3^ 
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many ways, etc.. 
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Hence we conclude 
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Putting these two together we find 
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This result has two conjugate poles in the complex density plane at a distance 



(23) 



C2 



P' <p, 



-(S)2Ar 

from the origin. It is clear that the model- and approximation-dependent sign of (S)2Ar 
determines the position of the poles in the complex p plane, and in particular if they are on 
the real axis, or not. For instance in the Born approximation to the linear sigma model 
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(24) 

Thus in the hnear sigina model (E)2Af is positive, as given in Eq. (|14D, and the poles are 
at imaginary densities ±ipc2 a factor \/2 larger in absolute value than pci predicted by the 
one-nucleon rescattering series. 

This sort of analysis can and ought to be extended to 3-, 4-, and higher N-pion rescat- 
tering graphs. We close with a conjecture: In the linear sigma model the 3-nucleon critical 
density pca is larger than the two-nucleon one pc2, the 4-N higher pc^ than the 3-N one, etc.. 

IV. SUMMARY AND CONCLUSIONS 

In summary, in this paper we have: (i) given an explicit formula for the n*^-order term 
in the powers-of-density (p") expansion of the nuclear matter sigma term (Zl)p and of the 
nuclear matter fermion condensate {qq)p in particular; (ii) calculated the 0{p^) coefficient 
in the Born approximation to the linear sigma model; (iii) iterated these primitive 0{p'^) 
terms to infinite order; (iv) found poles and radii of convergence associated with two such 
resummations. Thus we found that the "larger" of the two sums actually has the larger 
radius of convergence. 

Perhaps the most important conceptual contribution of this paper is that of putting the 
older "intuitive" calculations onto formally sounder grounds of QFT. In particular, we have 
shown that the fundamental elements of a calculation of (S)p are the connected n-nucleon 
sigma term matrix elements. All calculations of {qq) p can be reduced to these elements. 
Moreover, this result gives a clear definition of the "NN correlations" in this context. 

Some of the ideas introduced in this paper may have been tacitly assumed in earlier 



work, most notably that by the Manchester P] and the Lyon groups [Hllini- ^- Ericson 



correctly concluded that the resummation of the single-nucleon sigma term leads to a 
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geometric progression in the baryon density. Similarly, Birse and McGovern [Q came close 
to formulating the correct O(p^) prediction of the LDT in the linear sigma model. Finally, 
some of the diagrammatic methods used here were developed in the 60's and 70's P, p!3|JT^ 
for the description of pion propagation in nuclear matter. 
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FIG. 1. Feynman diagrams contributing to the "elementary" vrA^ elastic scattering amplitude: 
(a) the sigma-pole graph, and the (b), (c) the nucleon-pole graphs. The zig-zag line denotes a 
sigma meson, the dashed one a pion and the solid one a nucleon. 
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FIG. 2. Feynman diagrams contributing to the "elementary" tto" elastic scattering amplitude: 
(a) and the (b), ie., to the sigma-meson's sigma term. 
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FIG. 3. An example of order n = 3 reducible single-nucleon pion-rescattering graph. 
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FIG. 4. The pion-two-nucleon scattering amplitude: (a) the sigma-meson sigma term effective 
graph; (b) intermediate-state interaction graphs, (c) initial- and final-state interaction graphs. 
Graphs (b) and (c) may be termed NN correlation effects. The square box with four external 
nucleon lines denotes and NN potential due to the exchange of a single pion and sigma meson. 
The hatched "blob" in (a) represents the "elementary" irN elastic scattering amplitude, such as 
the one shown in Fig. (1) in the linear sigma model. 
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